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Abstract 

The generalized tight-binding model, with the exact diagonalization method, is 
developed to investigate optical properties of graphene in five kinds of external fields. 
The quite large Hamiltonian matrix is transferred into the band-like one by the re¬ 
arrangement of many basis functions; furthermore, the spatial distributions of wave 
functions on distinct sublattices are utilized to largely reduce the numerical compu¬ 
tation time. The external fields have a strong influence on the number, intensity, 
frequency and structure of absorption peaks, and the selection rules. The optical 
spectra in a uniform magnetic field exhibit plentiful symmetric absorption peaks and 
obey a specific selection rule. However, there are many asymmetric peaks and extra 
selection rules under the modulated electric field, the modulated magnetic field, the 
composite electric and magnetic fields, and the composite magnetic fields. 
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1 Introduction 


Monolayer graphene (MG), constructed from a single layer of carbon atoms densely packed 
in hexagonal lattice, was successfully produced by mechanical exfoliation. This partic¬ 
ular material offers an excellent system for studying two-dimensional (2D) physical prop¬ 
erties, such as the quantum Hall effects, |3|-[I^ and these properties could be preliminar¬ 
ily comprehended by the energy dispersion (or called energy band structure), which can 
directly reflect the main features of electronic properties. In the low-energy region of 
< 1 eV, MG possesses isotropic linear bands crossing at the K (K') point and is 
regarded as a 2D zero-gap semiconductor, where c (v) indicates the conduction (valence) 
bands. [Ib] The linear bands are symmetric about the Fermi level = 0) and become non¬ 


linear and anisotropic with > 1 eV. 16 Most importantly, the quasiparticles related 


to the linear bands can be described by a Dirac-like Hamiltonian, 17 which is associated 


with relativistic particles and dominates the low-energy physical properties. [^, 18,19 Such 
a special electronic structure has been verified by experimental measurements. [2, 20 

MG has become a potential candidate of nano-devices due to its exotic electronic prop¬ 
erties. Well understanding the behavior of MG under external fields is useful for improving 
the characteristics of graphene-based nano-devices. Five cases of external fields (see table 
below), which can be experimentally produced,are often applied to investigate the 
physical properties of few-layer graphenes (FLGs). In the presence of a uniform magnetic 
field (UM), the electronic states corresponding to the linear bands change into Landau lev¬ 
els (LLs) which obey a specific relationship E’^''^ oc where n'^ (rf) is the quantum 

number of the conduction (valence) states and Bum the magnetic field strength. The re¬ 
lated anomalous quantum Hall effects and particular optical excitations have been verified 


experimentally. 14,15 For a modulated magnetic field (MM), quasi-Landau levels (QLLs) 
possessing anisotropic behavior and the related optical absorption spectra with specified 


selection rules were shown. 28,29 Furthermore, Haldane predicted that MG in the modu- 
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lated magnetic field could reveal quantum Hall effects even without any net magnetic fiux 


through the whole space. 13 Concerning a modulated electric potential (ME), the linear 


dispersions become oscillatory and extra Dirac cones are induced by the potential. 30 -35 


Such a potential changes MG from a zero-gap semiconductor into a semimetal 33,35 and 


makes MG exhibit Klein paradox effect associated with the Dirac cones. For two cases 
of composite fields, a uniform magnetic field combined with a modulated magnetic field 
(UM-MM) and a uniform magnetic field combined with a modulated electric potential 
(UM-ME), the LL properties are drastically changed by the modulated fields. For both 
composite field cases, an unusual oscillation |36| - |40] of the density of states (DOS) similar 
to the Weiss oscillation obtained in 2D electron gas (2DEG) were shown. Furthermore, 
the broken symmetry, displacement of the center location, and alteration of the amplitude 


strength of the LL wave functions were also obtained. 41,42 


Graphene-related systems are predicted to exhibit rich optical absorption spectra. The 
spectral intensity of MG is proportional to the frequency, but no prominent peak exists at 


(u < 5 eV. 43 In FLGs, the interlayer atomic interactions drastically alter the two linear 


energy bands intersecting at Ep = 0. 44-54 As a result, conspicuous absorption peaks arise 


in optical spectra, 55 where the peak structure, intensity and frequency are dominated by 
the layer number and the stacking configuration. Furthermore, under an external perpen¬ 
dicular electric field or a uniform perpendicular magnetic field, the main features of the 
optical properties of the FLGs are strongly modified. |4^[5^[57| For theoretical studies, the 
complexity of calculating the optical absorption spectra is solved by the gradient approx¬ 
imation based on the generalized tight-binding model with exact-diagonalization method 
or effective-mass approximation. The way in which one can control the absorption peaks 
and selection rules is worthy to be reviewed in detail. 

On the other hand, there has been a considerable amount of experimental research on 
graphene-related systems under a uniform perpendicular magnetic field. From the mea¬ 
sured results, the features of MG and bilayer graphene are reflected in the magneto-optical 
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spectra. 58, 59 That is to say, the LL energies are proportional to or 

For any graphene system, the selection rule coming from the LLs close to Ep = 0 is 

An = = l.|60|-[6^ Moreover, similar results may also be found in AB-stacked 


graphite. 62 However, experimental measurements on optical properties under a non- 
uniform or composite helds are not available so far. 

In this chapter, we would like to focus on the optical absorption spectra of monolayer 
graphene under the hve cases of external helds, UM, MM, ME, UM-MM and UM-ME cases. 
The tight-binding (TB) model with exact diagonalization method is introduced to solve 
the energy dispersions and then the gradient approximation is applied to obtain the opti¬ 
cal absorption spectra. The main features of electronic properties, which include energy 
dispersions and wave functions, will be shown to comprehend the optical absorption prop¬ 
erties, where the dependence of absorption frequency on external helds, optical selection 
rules and anisotropic behavior will be discussed in detail. In Sec. 10.1, the tight-binding 
model corresponding to the hve cases of external helds is shown. In Secs. 10.2 to 10.6, the 
optical absorption spectra of MG under the UM, MM, ME, UM-MM and UM-ME cases 
will be reviewed, respectively. Einally, concluding remarks are presented in Sec. 10.7. 
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Physical properties of graphene under external fields 

External fields 

Related physical properties 

Uniform magnetic field 

Landau level and Abnormal quantum 

Hall effect, (fi -15 Magneto-optical selection 
rulep, 59 

Modulated magnetic field 

Quasi-Landau level, (2^ Quantum Hall 

effect without Landau level|13] 

Modulated electric potential 

Number increasement of Dirac cone,|30|- 

Semiconductor-metal transition, (M, 35 

Klein tunneling 17 

Uniform magnetic field+Modulated 

magnetic field, Uniform magnetic 

field+Modulated electric potential 

Weiss oscillation, 36 -40 Destruction of 

Landau-level wavefunction 41,42 
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2 Tight-Binding Model with Exact Diagonalization 


The low-frequency optical properties of graphene are determined by the vr-electronic struc¬ 
ture due to the 2pz orbitals of carbon atoms. The generalized tight-binding model with 
exact diagonalization method is developed to characterize the electronic properties and 
then the gradient approximation is applied to obtain the optical-absorption spectra. In the 
absence of external fields, there are two carbon atoms, the a and b atoms, in a primitive unit 
cell of MG, as shown in Fig. 10.1(a) by the green shadow, where the x- and y-direction are 
respectively the armchair and zigzag directions of MG. This indicates that the Bloch wave 
function T is a linear superposition of two TB functions associated with the 2pz orbitals 
and expressed as T = (pa ± (p;,, where ipa and (pb respectively stand for the tight-binding 
functions of the a and b atoms and are represented as 


= ^ exp(ik • Ra)x(r - R„), 

a 

( 10 . 1 a) 

(Pfe = exp(zk • Rfe)x(r - R^). 

( 10 . 1 b) 


x(r) is the normalized orbital 2pz wave function for an isolated atom. Moreover, the symbols 
7o (= 2.5 eV) and b' (= 1.42 A) shown in Fig. 10 . 1 (a) represent the nearest-neighbor 


atomic interaction (or called hopping integral) and the C-G bond length, respectively. 16 


Throughout this chapter, only 70 is taken into account and other atomic interactions are 
neglected. 

In the presence of an external field, the primitive unit cell is no longer the one shown in 
Fig. 10.1(a) since the external field leads to a new periodic condition. Here we choose the 
rectangular unit cell marked by the green rectangle in Fig. 10.1(b) as the primitive unit 
cell of graphene under the five kinds of external fields, where R = Rum, Rmm, Rme, and 
Rc (defined in the following) describe respectively the periods resulting from the uniform 
magnetic field, modulated magnetic field, modulated electric potential, and composite field. 
The major discussions are focused on R along the armchair direction. Consequently, an 
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Figure 10.1. The primitive unit cell of monolayer graphene (a) in the absence and (b) in 
the presence of external fields. 
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enlarged rectangular unit cell induced by an external field encompasses 2R a atoms and 2R 
b atoms. This implies that R determines the dimension of the Hamiltonian matrix, which 
is a 4i? X 4i? Hermitian matrix spanned by 4i? TB functions associated with the 2R a 
atoms and 2R b atoms. Based on the arrangement of odd and even atoms in the primitive 
cell, the Bloch wave function |Tk) can have the expression: 


2R-1 


2R 




( 10 . 2 ) 

m=l m=2 

kmk) (l^mk)) IS the TB function corresponding to the 2pz orbital of the mth a (b) atom. 
{A'^'") and are the subenvelope functions standing for the amplitudes of the wave 

functions of the a- and 6 -atoms respectively, where o (e) represents an odd (even) integer. 
Since the features of A^^'" {Bg'") and A^^'^ {B^'’) are similar, choosing only the amplitudes 
A^^'^ and B^’^’ is sufficient to comprehend the electronic and optical properties we would 
like to discuss in this chapter. The AR x 4R Hamiltonian matrix, which determines the 
magneto-electronic properties, is a giant Hermitian matrix for the external fields actually 
used in experiments. To make the calculations more efficient, the matrix is transformed into 
an M X AR band-like matrix by a suitable rearrangement of the tight-binding functions, 
where M is much smaller than AR. For example, one can arrange the basis functions 

as the sequence: |aik), \b 2 Ri,), Ihk), \a 2 Rk), |o 2 k}, | 62 ij-ik), | 62 k), |a 2 ii-ik}, .|oii-ik), 

\bR+ 2 k), \bR-ik), |oii+ 2 k}, \aRk), \bR+ik), \bRk); kii+ik)- Furthermore, distributions of the 
subenvelope functions are used to reduce the numerical computation time. The exact 
diagonalization method for numerical calculations is applicable to many kinds of magnetic, 
electric and composite fields. 


For the UM case Bum = Bum'z, a Peierls phase |28[|63 -65 related to the vector potential 
AuM=BuMxy is introduced in the TB functions. The phase difference between two lattice 
vectors (R^ and Rm') is defined as Gum = ^ Aum • dr , where 0o = hcje = 
4.1356 X 10 [T m^] is the flux quantum. The Peierls phase periodic along the armchair 
direction provides a specific period set as Rum = related Hamiltonian 
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is a ARum X ARum Hermitian matrix. The site energies, the diagonal matrix elements 
(amk|-f^|«mk) and (6mk|-f^|^mk), are set to zero and the nonzero matrix elements related to 
7o can be formulated as 

(^mk|.^^|nm'k) 7o®xpi[k • R-m') T Gum\- (10.3) 

Two kinds of periodic modulation fields along the armchair direction, the MM and 
ME cases, which can drastically change the physical properties of MG, are often se¬ 
lected for a study. For the MM case, Bmm = Rmm sin(27rx//MM) ^ is exerted on MG 
along the armchair direction, where Bmm is the field strength and 1mm is the period 
length with the modulation period Rmm = 1mm /^b'. The vector potential is chosen 
as Amm= {—Bmm^ cos{27rx/lMM))y and the corresponding Peirls phase is Gmm = 
Amm • dr. Thus the Hamiltonian matrix elements, which are similar to those in Eq. 
(10.3), are represented as 

7o®Xpi[k • (R,^ Rm') T GmmI- (10’^) 

For the ME case, Vme{x) = VMECOs(27rx/lME) along the armchair direction with the 
potential strength Vme and the period length Ime is taken into account. As the period is 
sufficiently large, the electric potential affects only the site energies but not the nearest- 
neighbor hopping integral. As a result, the site energies become 

H VMi?COs[(?7l X'Jtv/RmE^ — km; (10.5a) 

{bmk\H\bmli) = VMECOs[{m — 2/3)7 t/Rme] = Kn+1/3) (10.5b) 

where Rme = iME/^b' is the modulation period. The Hamiltonian matrices for the modu¬ 
lated magnetic field and the modulated electric potential are ARmm x ARmm and ARme x 
ARme Hermitian matrices, respectively. 

For a composite field case, a new periodicity, which is associated with periods induced 
by a uniform magnetic field and a modulated field, has to be defined. The rectangular unit 
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cell is enlarged along the x-direction and the dimensionality of the Hamiltonian matrix has 
to agree with the least common multiple of Rum and Rmm {Rum and Rme) for the UM- 
MM (UM-ME) case, namely Rq- The rectangular unit cell corresponding to each composite 
field contains 4i?e atoms {2Rc a atoms and 2Rc b atoms), and the magneto-electronic wave 
functions are linear combinations of the ARq TB functions. In a composite field case, the 
matrix elements are superposed by the elements associated with each combined external 
field. For the sake of convenience, we put the matrix elements in Eqs. (10.3)-(10.5) together 
as a common case and the elements are rewritten as 


(6^k|^^|om'k) = 7 oexpi[k ■ (R^-R^/) Gum + G'mm], (10.6a) 

(®mk|.^|®mk) ^mi (10.6b) 

{hmVi\H\bm\,) = Vm+\/Z- (10.6c) 


The off-diagonal elements are associated with the Peierls phases induced by the magnetic 
fields and the diagonal elements are related to the site energies induced by the modulated 
electric field. By diagonalizing the matrix, the energy dispersion E^’’^ and the wave function 
are obtained. It should be noted that the fc^j-dependent dispersions can be ignored 
when the period R is sufficiently large and thus only fc^-dependent dispersions are shown 
for the following discussions. 

When a monolayer graphene is excited from the occupied valence to the unoccupied 
conduction bands by an electromagnetic field, only inter-7r-band excitations exist at zero 
temperature. Based on the Fermi’s golden rule, the optical absorption function results in 
the following form 


A{uj) oc [ 

~ J \stBZ 


dk 


(4>^(k,n)|^|vl/"(k,n')} 


rUp 


X Im 


/(E^(k,n))-/(E’^(k,n')) 


(10.7) 


n) - E^{k, n') - uj - iTj ’ 
where /(£'(k, h)) is the Fermi-Dirac distribution function, and T (= 2 x 10“^7o) is the 
broadening parameter. The electric polarization E is the unit vector of an electric po- 
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larization. Results for E along the armchair and zigzag directions are taken into account 


for discussions. Within the gradient approximation, [66-68 the velocity matrix element 

h')) is formulated as 


2Rc 


^ (A" X B")Vi,(a„k|ff|6„.k) + h.c.. 


( 10 . 8 ) 


m,mf=l 


Equation (10.8) implies that the main features of the wave functions are major factors in 
determining the selection rules and the absorption rate of the optical excitations. Sim¬ 
ilar gradient approximations have been successfully applied to explain optical spectra of 


carbon-related systems, e.g., graphite, 69 graphite intercalation compounds, [70 carbon 


nanotubes, 71 few-layer graphenes, [^ and graphene nanoribbons. [22 
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3 Uniform Magnetic Field 


3.1 Landau Level Spectra 


In this section, we mainly focus on drastic changes of the Dirac cone as the result of a 
uniform perpendicular magnetic field. The magnetic field causes the states to congregate 
and induces dispersionless Landau levels, as shown in Fig. 10.2(a) for Bum = 5 T at 
fca, = 0. The unoccupied LLs and occupied LLs are symmetric about the Fermi level 
[Ep = ^). Each LL is characterized by the quantum number which corresponds to the 


the number of zeros in the eigenvectors of harmonic oscillator. 50,73 Each LL is fourfold 


degenerate without considering the spin degeneracy. Its energy may be approximated by a 
simple square-root relationship \E^'’\ oc ^h'^’^Bum, 


12 


74 


which is valid only for the range 


of < ±1 eV.^ 


3.2 Landau Level Wave Functions 

The LL wave functions, as shown in Eigs. 10.2(b) and 10.2(c), exhibit the versatility of 
spatial symmetry and can be described by the eigenvectors ((p„ (a:)) of harmonic oscillator, 
which obey the relationships, ((p„ {x) \ipn' {x)) = Sn,nf and ipn (a:^) = 0 for n < 0. The wave 
functions are distributed around the localization center, that is at the 5/6 position of the 
enlarged unit cell. Similar localization centers corresponding to the other degenerate states 
occur at the 1/6, 2/6, and 4/6 positions. The subenvelope functions can be expressed as 

(X (a:i) ± {X 2 ) , (^^l) T {X 2 ) , 

OC (Pnc^v_i {xs) ± (fnc^v {x^) , OC (^s) T ( 2 ^ 4 ) , 

for xi = 1/6, X 2 = 5/6, X 3 = 2/6, and X 4 = 4/6. (10.9) 

However, it is adequate to only consider any one center in evaluating the absorption spectra 
due to their identical optical responses. 


12 









0.05 


0.04 


0.03 


0.02 

^ 0.01 

o 

W 

0 

- 0.01 

- 0.02 


-0.03 

0 0.5 1 

ky (ttA/B b’) 


(a) n‘==3 
Yf=2 

rf=l 

Bua^5 T 

Ej^O 

n^=l 



0.81 0.82 0.81 0.82 
am{bm)/2R UM Clm (b„) / 2R UM 


Figure 10.2. The Landau level spectrum for the uniform magnetic field Bum = 5 T. The 
Landau level wave functions corresponding to (b) the a- and (c) the 6-atoms. 
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Figure 10.3. The optical absorption spectra for (a) Bum = 5 T and (b) Bum = 10 T. (c) 
The dependence of the absorption frequency on the square root of field strength Bum- 
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3.3 Optical Absorption Spectra of Landau Levels 


The low-frequency optical absorption spectrum of the LLs presents many interesting fea¬ 
tures as shown in Fig. 10.3(a) for Bum = 5 T. The spectrum exhibits many delta-function- 
like symmetric peaks with a uniform intensity. Such peaks suggest that LLs possess a 
zero-dimensional (OD) band structure or density of states. The optical transition channel 
with respect to each absorption peak can be clearly identified. A single peak is gener¬ 
ated by two transition channels n'LL—)• nLL and nLL—> n'LL, where the symbol n' —>■ n is 
used, for the sake of convenience, to represent the transition from the valence states with 
n' to the conduction states with n throughout this chapter. The quantum numbers related 
to the LL transitions must satisfy a specific selection rule, i.e.. An = \n'^ — n^| = 1. The 
selection rule is established by the main features of the wave functions. The velocity matrix 
M™, a dominant factor for the excitations of the prominent peaks, strongly depends on 
the number of zeros of and It has non-zero values only when and B^^\ 

expressed in orthogonality of ipn (x), possess the same number of zeros. Moreover, exam¬ 
ining all the transitions reveals the following relationship: oc B^'"{rf’'" + !)> with 

Aq = Aq and B^ — —B^. In other words, the quantum numbers of the conduction and 
valence LLs differ by one when Ao^^ and have the same ipn (3^)- 

In addition to the optical selection rules, the peak intensity and absorption frequency 
also deserve a discussion. In Fig. 10.3(b), the peak intensity is strengthened, whereas 
the peak number is reduced as the field strength increases. This is a result of the high 
degree of degeneracy in the first Brillouin zone and the expanded energy spacing between 
the LLs. The field-dependent absorption frequencies of the first four peaks a;|^, 

and are shown in Fig. 10.3(c). The frequencies become much higher in a stronger 
field. There exists a special square-root relation between co'2^ and Bq, i.e., co^u oc ^/Bum-, 
which has been confirmed by magneto-optical spectroscopy methods, such as experimental 
measurements of the absorption coefficient, mzn cyclotron resonance, |77|-[7^ and quantum 
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Hall conductivity.!^ 80-^ This square-root relation only exists in the lower frequency 
range lo < O.dqo (~1 eV). In the higher frequency range, LLs are too densely packed to be 


separated from one another. [I^ This leads to the disappearance of the relation between 
and Bum- 
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4 Spatially Modulated Magnetic Field 


4.1 Quasi-Landau Level Spectra 

Compared with the uniform case, a modulated magnetic held has a different impact on the 
electronic properties and leads to the diverse features observed in the optical absorption 
spectra. The presence of a modulated field has multiple effects on the energy bands, as 
shown in Fig. 10.4 for Bmm = 10 T and Rmm = 500. In the lower energy region, 
parabolic subbands appear around ky = ki = 2/3. The conduction and valence subbands 
are symmetric about the Fermi level {E^ = 0). The subbands nearest to Ep = 0 are 
partially hat and nondegenerate. The other parabolic subbands characterized by weak 
energy dispersions have double degeneracy and one original band-edge state at k^. The 
modulation effects on parabolic energy subbands result in four extra band-edge states at the 
sites on both sides of ki. They demonstrate the strongest dispersion and destruction of the 
double degeneracy. The low-energy subbands are regarded as quasi-Landau levels, which 
exhibit similar features of the LLs generated from a uniform magnetic held. Moreover, the 
ky range with respect to the weak dispersion and partial hat bands grows with increasing 
held strength and a longer modulation period. On the contrary, when the inhuence of the 
modulation held become much weaker with increasing energy, the parabolic subbands in 
the higher energy region are similar to the twofold degenerate subbands directly obtained 
from the zone folding of MG in the Bmm = 0 case (not shown). 

4.2 Quasi-Landau Level Wave Functions 

In the presence of a modulated magnetic held, the alterations of the wave functions are 
rather drastic. First, the QLL wave functions corresponding to ki are shown in Figs. 
10.5(a)-(f). The wave functions are composed of two tight-binding functions centered at 
Xi and X 2 . (B^) has two subenvelope functions A'^{xi) {B^{xi)) and A^(x 2 ) {B^{x 2 )) 
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Figure 10.4. The energy dispersions and the illustration of optical excitation channels for 
the modulated magnetic field along the armchair direction with Rmm = 500 and Bmm = 10 

T. 
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centered at = 1/4 and 0:2 = 3/4 of the primitive unit cell, respectively. The positions 
Xi and X 2 are located at where the field strength is at a maximum. The number of zeros 
of A^(x 2 ) {B^{xi)) is higher than that of A^{xi) {B^{x 2 )) by one at each QLL. A similar 
behavior is also shown by the valence wave function, where only the sign is flipped in either 
Aq or B^. The effective quantum number is defined by the larger number of zeros of 
the subenvelope functions. In addition, the twofold degenerate QLLs have similar wave 
functions (black curves and red dashed curves), with the only difference in terms of the 
sign change in the subenvelope functions. The wave functions at ki can be expressed as 

^ (X (xi) ± {X 2 ) , a;, (X (a^i) t (^ 2 ) , 

g oc (xi) T ( 2 ^ 2 ), -B/e OC (xi) ± (X 2 ), 

for xi = 1/4 and X 2 = 3/4. (10.10) 

The wave functions would be strongly modified as the wave vectors gradually move away 
from ki. Secondly, the wave functions at several special k points are illustrated to examine 
the effects caused by the modulated magnetic field. As the wave vector moves to ky = k 2 , 
the doubly degenerate QLL starts to separate into two subbands. The two subenvelope 
functions A^{xi) (B^{xi)) and A^{x 2 ) {B^{x 2 )) move toward each other and shift to the 
center of the primitive unit cell with nearly overlapping, as shown in Fig. 10.5(g) and 
10 .5(h). At ks and k^, the higher and lower subbands have the extra band-edge states la 
and ip, respectively. The subenvelope functions of the la state, as shown in Fig. 10.5(i) 
and 10.5(j), exhibit a strong overlapping behavior compared to those at ky = k 2 (red-dashed 
curves in Fig. 10.5(g) and 10.5(h)). Similar behavior can also be found in the wave functions 
at 1/5. This implies that there is a higher degree of overlap in the subenvelope functions at 
the extra band-edge states rf''^a and n‘^’^/5. Moreover, the two states associated with the 
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different linear combinations of A^(xi) {B^{xi)) and A^{x 2 ) {B^{x 2 )) are represented as 


^o,e ^n=-l (^l) + (^2) 

B^^ OC (xi) - {X2) 

OC d'n.-i (xi) - {X2) 

-®o,e ^n” (3^1) + ^n”-l (^2) 

for Xi ~ X 2 — 1/2. 


for 72*^0; and (xi) — d'^c (X2) for n‘^/5, 

for 72*^0; and \l/„c (xi) + ^nc_i (X 2 ) for rfl3, 
for rAa and ^nf-i (xi) + \l/„« (X 2 ) for rA(3, 
for rAa and d'„t> (xi) — (^^ 2 ) for rfP, 


( 10 . 11 ) 


4.3 Optical Absorption Spectra of Quasi-Landau Levels 

Under the modulated magnetic field, the parabolic energy bands possess several band- 
edge states. A wave function composed of two tight-binding functions presents a complex 
overlapping behavior. The above-mentioned main features of the electronic properties 
are expected to be directly reflected in optical excitations. The low-frequency optical 
absorption spectra for Rmm = 500 and Bmm = 10 T, as shown in Fig. 10.6(a) by the black 
and blue solid curves for E _L x and E || x respectively, exhibit rich asymmetric peaks in the 
square-root divergent form. These peaks can be divided into the principal peaks cup’s and 
the subpeaks cu^’s based on the optical excitations resulting from the original band-edge 
and extra band-edge states, respectively, cu^’s can be further classified into two subgroups 
Ug’s and a;|’s which primarily come from the excitations of extra band-edge states a ^ {3 
{(3 —7- a) and a ^ a (P ^ P), respectively. What is worth mentioning is that the 
spectra for E T x and E || x are distinct, especially for the subpeaks cu^’s. The former is 
mainly composed of the subgroup cug, while the latter mainly consists of the subgroup a;|. 
This implies that the optical absorption spectra reflect the anisotropy of the polarization 
direction. For the modulation along the zigzag direction at Rmm = 866 and Bmm = 10 T, 
the absorption spectrum (red dashed curve in Fig. 10.6(a)) shows features similar to those 
of the spectrum corresponding to the armchair direction at Rmm = 500 and Bmm = 10 
T. Rmm = 866 for the zigzag direction and Rmm = 500 for the armchair direction possess 
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Clm {bm) ! '^RmM 


Figure 10.5. The wave functions of Qnasi-Landan levels at (a)-(f) the original band-edge 
state ki with the quantnm numbers rf = 1 and rf = 2, (g) and (h) the split point 

^2 with = 1 , and (i)-(l) two extra band-edge states and ^4 with = 1 . 
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the same period length based on the definitions Rmm = 1mm /^b' and Rmm = 1mm /V^b' 
associated with the zigzag and armchair directions, respectively. Moreover, the anisotropic 
features of the modulation directions will be revealed in the higher frequency region or the 
smaller modulation length. 

As the field strength rises, the peak height and frequency of the principal peaks increase, 
and the peak number decreases, as shown in Fig. 10.6(b) by the red curve for Rmm = 500 
and Bmm = 20 T. These results mean that the congregation of electronic states is more 
pronounced as the field strength grows. In addition to the field strength, the optical- 
absorption spectrum is also influenced by the modulation period. In Fig. 10.6(b), the blue 
curve shows the optical spectra of Bmm = 10 T for Rmm = 1000. The subpeaks strongly 
depend on the period, i.e., they represent different peak heights and frequencies with the 
variation of Rmm- However, the opposite is true for the principal peaks. 

The peaks in the low-frequency absorption spectra can arise from the different selection 
rules. Fig. 10.4 illustrates the transition channels of the principal peaks resulting from the 
original band-edge states denoted as cup^’s in Fig. 10.6(a). Each ujpn corresponds to the 
transition channels from QLLs n —)• n -|- 1 and n —)• n -|- 1 at the original band-edge state 
and the selection rule is represented by An = \rf' — 'n}’\ = 1 which is same as that related 
to LLs. The main reason for this is that the subenvelope functions {A'^^\x 2 )) 

and Bo^‘^\xi) {Bo^‘^\x 2 )) associated with the effective quantum numbers n -|- 1 (n) and n 
(n + 1) have the same number of zeros, respectively. As discussed in the former section, 
peaks arise in the optical absorption spectra when the number of zeros is the same for Ao 
and Bo^^^ in Eq. (10.10). The subpeaks originating from the extra band-edge states display 
a more complex behavior. The excitation channels for the subpeaks and in Eig. 
10 .6(a) are shown in Eig. 10.4. The subpeaks of different selection rules. An = 0 and 1, 
come into existence simultaneously. Eor example, 00^2 comes from the excitation channel 
la —7- ip (1/3 ^ la) and comes from the excitation channel 1/3 ^ 2a (2a —?• 1/3). The 
extra selection rule An = 0 reflects the overlap of subenvelope functions A'^{xi) (B^{xi)) 
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Figure 10.6. The optical absorption spectra for (a) Bmm = 10 T at a fixed periodic 
length with modulation and polarization along the armchair and zigzag directions and (b) 
different modulation periods and field strengths with both the modulation and polarization 
along the armchair direction. 
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and A^(x 2 ) {B^{x 2 )) located around a:i ~ X 2 ~ 1/2. The subenvelope functions A'^^\xi) 
{A'^^\x 2 )) and BZ^^\x 2 ) {B'^^\x\)) of the effective quantum number n also have the same 
number of zeros at the identical position, a cause leading to the extra selection rule An = 0. 

The frequency of principal peaks in the optical absorption spectra is worth a closer 
investigation. The relation between the frequencies of the first four principal peaks and the 
modulation period is shown in Fig. 10.7(a). The wp’s present a very weak dependence on 
the period as Rmm becomes sufficiently large, whereas they exhibit a strong dependence on 
the field strength. The frequencies grow with increased Bmm, as shown in Fig. 10.7(b). The 
dependence of (Up’s on Bmm is similar to what is seen in the case of a uniform perpendicular 
magnetic field, i.e., wp’s oc \/Bmm, as indicated by the red lines. The predicted results 


could be verified by optical spectroscopy. (14, 20,79 
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Figure 10.7. The dependence of the absorption frequency on (a) the period Rmm and (b) 
the square root of field strength Bmm- 
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5 Spatially Modulated Electric Potential 

5.1 Oscillation Energy Subbands 

Besides the spatially modulated magnetic field, the low-energy physical properties can also 
be strongly tuned by a spatially modulated electric potential. The energy bands for Vme = 
0.05 7o and Rme = 500 are shown in Fig. 10.8. The unoccupied conduction subbands 
are symmetric to the occupied valence subbands about Ep. The parabolic subbands are 
nondegenerate and oscillate near ky = 2/3. There exists on intersection where two parabolic 
subbands cross each other at Ep. Each subband has several band-edge states, which lead 
to the prominent peaks in the DOS and optical absorption spectra. For convenience, 
these band-edge states are further divided into two categories called fi and u states, as 
indicated in Fig. 10.8. The two ji (z/) states at the left- and right-hand sites of ky = 2/3 
might have a small difference in energies; that is, parabolic bands might be bilaterally 
asymmetric about ky = 2/3. Not far away from ky = 2/3, the energy subbands with linear 
dispersions intersect at Ep, preserving more Fermi-momentum states and forming several 
Dirac cones. Moreover, the number of Fermi-momentum states or Dirac cones increases 
with the potential strength and modulation period. 

The optical absorption spectrum for Rme = 500 and Vme = 0.05 70 along the armchair 
direction, as shown in Fig. 10.9 by the black solid curve, exhibits two groups of prominent 
peaks, Eji’s and T„’s. They are mainly due to the optical excitations from /i^ to lE^+i 
to and to iEn +2 {t^n +2 to /U^), respectively. Moreover, with regard to the peak 
intensity, the peaks E„’s (T„’s) can be further divided into two subgroups. For example, 
the peak heights of Ei, E3; E5, respectively, resulting from the transitions of ii\ to /i^ 
to /r^), /U3 to (iJ,^ to /Tg); nl to /Ug (/Ug to /Xg) are very low, while the peaks E2, E4; Eg 
originating from the excitations to /ig (/ig to /r^), jEl to /ig (/ig to /r^); /ig to /i^ (/iy to 
/ig) present much stronger intensities than the peaks Ei, Eg and Eg. That is to say, the 
peak of E2n’s are higher than those of E2n-i’s in the group E„. The peaks of T„’s exhibit 
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Figure 10.8. The energy dispersions for the modnlated electric potential along the arm¬ 
chair direction with Rmm = 500 and Vme = O.OSyo. 
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similar features to those of E^’s. For instance, peaks Ti, T3; T5, respectively arising from 
the transitions of to /Ug (/Ug to /r^), /Ug to fxl {nl to /Xg); nl to (/r^ to /U 5 ), own the 
peaks with very weak intensities. In contrast to Ti, Tg; T5, the peak intensities of T 2 and 
T 4 resulting from the excitations to (nl to /r^) and nl to /Xg (/Ug to /r^) are relatively 
stronger. Furthermore, the ji and z/ states lead to different contributions to the two kinds 
of optical absorption peaks. Most peaks originating from the two different band-edge states 
have nearly the same frequencies, while the peak intensities are not the same. The blue and 
red curves correspond to the optical absorption spectra which contains only the excitations 
of /i and u states, respectively. Except for the peak E 2 with comparable contributions 
which are attributed to the transitions of ji and f states, the other peaks with different 
contributions from the two states have nearly the same frequency. The peaks from the /r 
states exhibit much stronger intensities than those from the v states. In other words, peaks 
in the optical absorption spectrum mainly result from excitations of the ji states. 

5.2 Anisotropic Optical Absorption Spectra 

The polarization direction and the strength, period and direction of the modulating electric 
field strongly affect the features of the optical absorption spectrum. The spectra associated 
with E T T (black solid curve) and E || T (red solid curve) for Rme = 500 and Vme = 0.05 
7 o along the armchair direction and Rme = 866 and Vme = 0.05 70 along the zigzag 
direction (blue solid curve) are shown in Fig. 10.10(a) for a comparison. Compared with 
the results of E T T and E || T, the peak structures related to the two polarization 
directions are totally different, which reflect the anisotropic behavior of the polarization 
direction. Similarly, the anisotropy of the modulation directions is reflected by that the 
absorption spectra corresponding to the armchair and zigzag directions display distinct 
features, i.e., the anisotropic behavior of the polarization directions are more obvious than 
that in the MM case. With increasing the modulation strength to Vme = O.I 70 (red solid 
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Figure 10.9. The optical absorption spectra corresponding to Fig. 10.8, which includes 
the contributions from the /r and u states, respectively. 
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curve in Fig. 10.10(b)), the results show that the peak intensity strongly depends on Vmei 
but their relationship is not straight forward. For the modulation period, the spectra at a 
larger Rme = 1000 (blue solid curve) along the armchair direction present features diverse 
to those in the spectra at Rme = 1000 . The peak number grows and the peak intensities 
decay with a increase of the period. A redshift occurs in longer periods. For example, the 
peak frequencies Ei, E 3 ; E 5 , as indicated in black and green curves, are almost reduced to 
half of the original ones when the modulation period is enlarged from 500 to 1000. 

The optical absorption spectra in the ME case do not reveal certain selection rules. 
This is due to the fact that the amplitudes and of the wave functions do not exist 
a simple relationship similar to that in the UM and MM cases. The wave functions in the 
modulated electric potential are no longer distributed around the center location; rather, 
they display standing-wave-like features in the primitive unit cell and are distributed over 
the entire primitive cell, as shown in Fig. 10.11. However, the wave functions of the 
edge-states ji and e exhibit irregular behavior such as disordered numbers of zero points, 
asymmetric spatial distributions, and random oscillations. These irregular waveforms might 
result from different site energies for the carbon atoms in the modulated electric potential. 
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Figure 10.10. The optical absorption for (a) Vme = O.O 570 at a fixed periodic length with 
modulation and polarization along the armchair and zigzag directions and (b) different 
modulation periods and field strengths with both the modulation and polarization along 
the armchair direction. 


31 























Figure 10.11. The wave functions at different band-edge states, the /U? and I'f states for 
i = 1'7. 
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6 Uniform Magnetic Field Combined with Modulated 
Magnetic Field 

6.1 Landau Level Spectra Broken by Modulated Magnetic Fields 

A further discussion of graphene in a composite field, the UM-MM case, is presented in this 
section. The main characteristics of the LLs at Bum = 5 T are affected by the modulated 
magnetic field {Bmm = 1 T and Rmm = 395), as shown in Fig. 10.12(a) by the black 
curves. The LL with = 0 at Ep = 0 remains the same features of the UM case. On the 
other hand, each dispersionless LL with ^ 1 splits into two parabolic subbands with 
double degeneracy. The subbands possess two kinds of band-edge states, and 
which correspond to the minimum field strength Bum — Bmm and maximum field strength 
Bum + Bmm, respectively. The surrounding electronic states at congregate more 

easily, which results in the smaller band curvature. Comparably fewer states congregate at 
and the resulting band curvature is larger. Increasing Bmm induces more complex 
energy spectra, as shown in Fig. 10.12(b) for Rmm = 395 and Bmm = 5 T. The parabolic 
subbands with ^ 1 display wider oscillation amplitudes, stronger energy dispersions, 
and greater band curvatures. The largest and smallest band curvatures occur at the local 
minima and local maxima states, respectively. The subband amplitudes are 

nearly linearly magnified by Bmm as Bmm < Bum- It is noticeable that neither the 
minima of the conduction bands nor the maxima of the valence bands exceed Ep = 0 even 
for Bmm much larger than Bum, as shown in Fig. 10.12(c) for Rmm = 395 and Bmm = 40 
T. Thus no overlap exists between the conduction and valence bands, regardless of the 
modulation strength. With further increasing modulated field strength as Bmm ^ Bum, 
the electronic structures are expected to approach to those in the MM case. 
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Figure 10.12. The energy dispersions for (a) the nniform magnetic field Bum = 5 T by the 
red curves and the composite field Bum = 5 T combined with Rmm = 395 and Bmm = 1 
T by the black curves, (b) Bum = 5 T combined with Rmm = 395 and Bmm = 5 T, and 
(c) Bum = 5 T combined with Rmm = 395 and Bmm = 40 T. All modulated fields are 
applied along the armchair direction. 
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6.2 Symmetry Broken of Landau Level Wave Functions 

The LL wave functions modified by the modulated magnetic held are shown in Fig. 10.13. 
The spatial distributions corresponding to labeled in Fig. 10.12, exhibit slightly broad¬ 
ened and reduced amplitudes, as indicated by the black curves in Figs. 10.13(a)-(d) for 
Bmm = 1 T. However, the spatial symmetry and the location centers of the wave func¬ 
tions remain unchanged. Under the inhuence of a small Bmm, the simple relation between 
and of the wave functions is almost preserved. However, a stronger modula¬ 
tion strength results in greater spatial changes of the wave functions, as shown in Figs. 
10.13(e)-(f) for Bum = Bmm = 5 T. The increased broadening and asymmetry of the 
spatial distributions of the wave functions at = 0 are revealed. However, the spatial 
distributions with ^ 1 are only widened (i.e., = 1 in Figs. 10.12(g) and (h)), but 

the spatial symmetry is retained. With increasing Bmm, as shown in Figs 10.13(i)-(l) for 
Bmm = 40 T, the symmetry of the wave functions with = 0 is recovered and one can 
expect that the main features of the wave functions will become similar to those in the 
MM case. Obviously, the electronic properties show critical changes as Bmm equals Bum, 
which should be reflected to the optical properties. 

6.3 Magneto-Optical Absorption Spectra with Extra Selection 
Rules 

The optical absorption spectra corresponding to Fig. 10.13 are shown in Figs. 10.14 and 
10.15. In Fig. 10.14(a), the absorption spectra corresponding to the UM-MM case at 
Bum = 5 T with Rmm = 395 and Bmm = 1 T and the UM case at Bum = 5 T are shown 
together for a comparison. The red curves coming from the LLs at Bum = 5 T display 
delta-function-like peaks uj^u with the selection rule An = 1. However, the modulated 
magnetic field modihes each delta-function-like peak into two split square-root-divergent 
peaks, and as shown by the black curves. Each originates from the 
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Figure 10.13. The wave functions with = 0 and = 1 at for (a)-(d) the uniform 
magnetic field Bum = 5 T by the red curves and the composite field Bum = 5 T combined 
with Rmm = 395 and Bmm = 1 T by the black curves, (e)-(h) Bum = 5 T by the red 
curves and Bum = 5 T combined with Rmm = 395 and Bmm = 5 T by the black curves, 
and (i)-(l) Bum = 5 T combined with Rmm = 395 and Bmm = 40 T. 
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transitions of —?• n + and n + —?• (nr] n + Ir] and n + It] ^ nrj) and 

its absorption frequency is same as that generated from the LLs at Bum — Bmm = 4 T 
(Bum — Bmm = 6 T). These absorption peaks obey a selection rule, An = 1, similar to 
that in the UM case. 

With increasing the modulated field strength as Bmm = Bum = 5 T, the absorption 
spectrum has evident variety, as shown in Fig. 10.14(b). In addition to the peaks n;””' and 
a;””' with the selection rule An = 1, two extra peaks with An = 2 and 3, and 
are generated. These two peaks do reflect the fact that the wave functions of the LLs with 
n*^’^ = 0 are destroyed by the modulated magnetic field. As the modulated field strength 
further raises to Bmm = 40 T (red dashed curves Fig. 10.15), the spectrum displays some 
features similar to those of the spectrum in the MM case at Bmm = 395 and Bmm = 40 T 
(black solid curves in Fig. 10.15), i.e., the principal peaks up’s and the subpeaks us’s in 
the MM case are also shown in the UM-MM case as Bmm > Bum- Moreover, the subpeaks 
features, which are associated with the positions at the net field strength equal to zero, 
are almost the same in both the MM and UM-MM case. The principal peaks, however, 
possess a pair structure with Up^ and which respectively correspond to two different 
field strengths, \Bum — Bmm\ = 35 T and \Bum + Bmm\ = 45 T, and thus the difference 
between two field strengths lead to distinct absorption frequencies. For Bmm ^ Bum, one 
can anticipate that the frequency discrepancy between the pair cop^ and Up^ becomes very 
small and then they will merge into one single peak, u)pn, i.e., the absorption spectrum 
restores to that in the pure MM case. 

The dependence of the absorption frequency on the modulated field strength is shown 
in Fig. 10.16 for Bmm < 5 T. In the range of Bmm A Bum, each of absorption peaks 
and n;””' is linearly dependent on Bmm- This reflects the fact that the subband amplitudes 
are nearly linearly magnified by Bmm within the range. However, in the higher absorption 
frequency region or the field range of Bmm > Bum, the linear-dependence relationship 
will be broken since the subbands become overlapping and the subband amplitudes are not 
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Figure 10.14. The optical absorption spectra corresponding to (a) the uniform magnetic 
field Bum = 5 T by the red curve and the composite field Bum = 5 T combined with 
Rmm = 395 and Bmm = 1 T by the black curve and (b) the composite field Bum = 5 T 
MM — 395 and Bmm = 5 T. 
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Figure 10.15. The optical absorption spectra corresponding to the composite held Bum = 
5 T combined with Rmm = 395 and Bmm = 40 T by the red dotted curve and the pure 
modulated magnetic held Rmm = 395 and Bmm = 40 T by the black curve. 
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linearly magnified by Bmm anymore. 
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7 Uniform Magnetic Field Combined with Modulated 
Electric Potential 

7.1 Landau Level Spectra Broken by Modulated Electric Poten¬ 
tials 

Compared with the situation in a modulated magnetic field, a modulated electric potential 
creates distinct effects on the LLs, as shown in Fig. 10.17(a) by the black curves for . The 
OD LLs become the ID sinusoidal energy subbands when electronic states are affected by 
the periodic electric potential. Each LL with four-fold degeneracy is split into two doubly 
degenerate Landau subbands (LSs), as shown in Fig. 10.17(a) by the black curves. Each 
LS owns two types of extra band-edge states, and For the conduction (valence) 
LSs, the band-edge states of and {k^^ and /c^) are, respectively, related to the 
wave functions, which possess a localization center at the minimum and maximum electric 
potentials (discussed in the ME case). It should be noted that the two LSs of rf’’" = 0 only 
have the p-type band-edge states. Under a small modulation strength, the energy spacings 
{Egs) or band curvatures for both k^^ and k^^ are almost the same, where Eg is the spacing 
between a LS and a LL at k^^ or k^^. On the other hand, when the modulation strength 
is sufficiently large (e.g., Vme = 0.02 70 in Fig. 10.17(b)), the energy dispersions of LSs 
are relatively strong and Eg decreases with increasing state energies, i.e., the oscillations 
of LSs decline with an increase of In comparison to the k^^ state, the state owns 
the smaller energy spacing and band curvature. Such differences are associated with the 
localization of the wave function within the potential well. 
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Figure 10.17. The energy dispersions for (a) the uniform magnetic field Bum = 5 T 
by the red curves and the composite field Bum = 5 T combined with Rmm = 395 and 
Vme = 0.00570 by the black curves, (b) Bum = 5 T combined with Rmm = 395 and 
Vme = 0.0270. All modulated fields are applied along the armchair direction. 
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7.2 Landau Level Wave Functions Broken by Modulated Electric 
Potentials 

The main features of the LL wave functions are altered by the modulated electric potential. 
The illustrated wave functions at the band-edge states and the midpoint {kmid, indicated in 
Fig. 10.17(a)) between two band-edge states are used to examine the modulation effects. 
The wave functions at kmid are modified by Vme{x), as shown in Figs. 10.18(a)-(h) for 
Bum = 5 T and Rme = 395 at Vme = 0 , 0.005 and 0.02 70 . of n'^'^ = 0 is slightly 
reduced, while B^’’" of = 0 is slightly increased (Figs. 10.18(a) and 10.18(b)) after Vme 
is introduced. This means that carriers are transferred between the a- and 6 -sublattices. 
With an increasing the spatial distribution symmetry of the LL wave functions is 
broken. The conduction and valence wave functions are, respectively, shifted toward the 
-[-X and —X directions, as shown in Figs. 10.18(e)-(h) for example. The proportionality 
relationship between of and of -|- 1 no longer exists, and neither do the 
relationships and B^ = —B^. Moreover, the stronger Vme leads to greater 

changes in the spatial distributions of the wave functions. The spatial distributions of the 
wave functions strongly depend on ky. As for the band-edge states, the aforementioned 
relationships of the wave functions are absent when n'^’^’s are sufficiently large enough. 
For small (Figs. 10.18(i)-(l)), wave functions are less influenced by the modulated 

electric potential. However, the spatial distribution of LS with a larger becomes wider 
(narrower) for the k^^ (fc^) state, as shown in Figs. 10.18(m)-(p). Moreover, the localization 
center of the band-edge states is hardly affected by Vme{x). 

7.3 Magneto-Optical Absorption Spectra Destroyed by Modu¬ 
lated Electric Potentials 

Under a modulated electric potential, the changes in the electronic properties of LLs are 
manifested in the optical absorption spectra. Each LL is split into two kinds of sinusoidal 
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Figure 10.18. The wave functions for (a)-(d) = 0, = 1, = 5, and = 5 at 

kmid, (i)-(n) = 0, = 1 and = 5 at and (l)-(m) = 5 at The results 

corresponding to the uniform magnetic field Bum = 5 T, the composite field Bum = 5 
T combined with Rmm = 395 and Vme = O.OO 570 , and Bum = 5 T combined with 
Rmm = 395 and Vme = O.O 270 are indicated by the red, black, and blue curves respectively. 
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subbands, with one leading the other by 1/6 of a period (Fig. 10.17(a)). The spatial local¬ 
ization region of the wave function is completely different between two kinds of subbands 
(not shown), but identical in corresponding to the same kind of LSs with different quantum 
numbers at the same ky. This indicates that the optical transition between two different 
kinds of subbands is forbidden. The absorption spectrum for Rme = 395 at Vme = O.OOSyo 
is shown in Fig. 10.19(a) by the black line. Each absorption peak, cu””/ originates from a 
transition between two LSs with quantum numbers n and n' for the same kind of subbands. 
In addition to the original peaks, which correspond to the selection rule An = 1 similar 
to that of LLs, there are extra peaks not characterized by the same selection rule. In the 
frequency range u < 0.1 70 , the peak intensity of significantly reduced by Vme, 

declines as the frequency increases. The extra peaks with An 7 ^ 1 behave the opposite way. 
For a small Vme, the peaks with An = 1 are much stronger than those with An 7 ^ 1 . 

The original and extra peaks can be explained by the subband transitions associated 
with a certain set of ky points. For the k^^^^ —>■ k^^^^ transitions, the corresponding band- 
edge states have a high DOS and the symmetry of wave function is little changed. Therefore, 
they can promote the prominent peaks consistent with the An = 1 selection rule. As shown 
in Fig. 10.19(b) by the thin dashed lines, and represent the absorption peaks 

from the transition channels [nfc^ ^(n -|- l)k^^, (n + i)k§^ “t(n -|- l)fc^, 

(n -|- 1)A:^ —>■ nk^^], respectively. These two peaks are close to each other, and almost 
overlap with the original peaks. But in cases where either n or n' is zero, only the peak 
can be created by two transition channels: [Ofc^g ^ Ifc^, Ifc^ ^ Ofc^g]. The reduction 
of the transition channels is due to the fact that the = 0 subbands oscillate between 
the conduction and valence bands. For the k^id kmu transitions, the significant change 
to the symmetry of the wave function results in different selection rules, i.e., An 7 ^ 1 . As 
shown in Fig. 10.19(b) by the thick dashed lines, the absorption peak corresponds 
to the transitions between two LSs of and from the k^id states. These types of 
peaks are responsible for the extra peaks in Fig. 10.19(a), i.e., the An 7 ^ 1 peaks primarily 
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Figure 10.19. The optical absorption spectra corresponding to (a) the uniform magnetic 
field Bum = 5 T by the red curve and the composite field Bum = 5 T combined with 
Rmm = 395 and Vme = O.OOSyo by the black curve, (b) the transitions of 
and kmid kmid for Bum = 5 T combined with Rmm = 395 and Vme = O.OOSyo by the 
dotted and dashed curves respectively, and (c) the higher frequency region of 0 . 1 ~ 0 . 27 o. 
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arise from the middle states at lower modulation strength. The k^id —?• kmid transitions also 
contribute to the original peaks. Such contributions decline with an increased frequency. 
However, the An = 1 absorption peaks are dominated by both the band-edge and middle 
states. 

The peak intensities from different selection rules change dramatically with respect 
to the variation in frequency. The intensities of the An = 1 peaks gradually rise for a 
further increase in frequency, while the opposite is true for those of the An ^ 1 peaks (Fig. 
10.19(c)). Within the frequency range of 0.1 70 < w < 0.2 70 , the former is lower than 
the latter. The main reason for this is the fact that the symmetry violation of the wave 
functions is enhanced for LSs with larger n'^’^^’s. It is deduced that the extra absorption 
peaks of An 7 ^ 1 are relatively easily observed for experimental measurements at higher 
frequencies. Each of them is composed of two peaks, which satisfy 

An = 0 and An = 2 with nearly the same frequency. 

The absorption spectrum exhibits more features for a higher modulated potential, as 
shown in Fig. 10.20 for Rme = 395 at Vme = 0, O.O 270 . In Fig. 10 . 20 (a), the intensity 
of the extra peaks becomes comparable to that of the original peaks. Unlike in the lower 
Vme case (Fig. 10.19(a)), the peak intensities, regardless of their types, vary irregularly 
with the frequency. The oscillations of the subband structure are obviously augmented 
and cause the subband transitions to change greatly with respect to ky. As a result, the 
absorption peaks grow much wider and split more evenly. Peaks generated by different 
selection rules are likely to appear, owing to the severe breakdown of the spatial symmetry 
of the wave functions. It should be noted that it is difficult to distinguish the original from 
the extra peaks solely based on the peak heights. Besides the aforementioned spectrum 
analysis, further discussions are made regarding two specific ky points. For the 
transitions, x- and p-type band-edge states show different behavior in terms of energy 
spacing and curvature such that the two transition channels -^{n + l)k^^, {n + l)k^^ 

and [nk^^ -^(n -\- 1 )A:^, (n -l- 1 )^^ ^ possess distinct frequencies (thin dashed 
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line in Fig. 10.20(b)). Therefore, the initially coinciding peaks split up, and the original 
intensities are divided into fractions. Moreover, the band-edge states can induce the extra 
peaks of An = 0, e.g., and When ky = kmid, fhe symmetry of the wave functions 
is destroyed. Consequently, the extra peaks are strengthened, and the original peaks are 
weakened or even disappear (thick dashed line in Fig. 10.20(b)). The corruption of the 
orthogonality of the sublattices and enables the subband transitions to occur 

from An = 3; such examples are seen in and The very strong dispersions of LSs 

also lead to the splitting of the original peaks associated with kmu- This can account for 
the ufj peaks and the ujH) peaks in Fig. 10.20(a). 

The modulation period strongly affects the magneto-optical spectrum of the LSs, while 
the modulation effect is less signihcant with sufficiently larger period. The absorption 
peaks from An = 1 are much higher than those from An 1 , as shown by the green 
line in Fig. 10.20(c) for Vme = 0.02 70 at a larger period R = 1580. The transition 
energies between the valence and conduction LSs at different ky values are almost the 
same. Thus, the frequencies associated with these LS tran¬ 

sitions are all the same. Since the symmetry of the wave function does not degrade 
much, the absorption peaks of An = 3 no longer exist. The modulated electric held 
E = —VxVMEi'x) = { 27 rVME/^b'RME)sm{ 27 rx/ 3 b'RME)x implies that the same value of 
Vme/Rme produces the same modulation effect. For instance, the absorption spectrum for 
Rme = 395 at Vme = O.OOSyo (the purple dashed curve in Fig. 10 . 20 (c)) is almost same as 
that for Rme = 395 at Vme = 0 , O.OOSyo owing to Vme/Rme = 0.005/395 = 0.02/1580. 

We look at the relationship between the absorption frequency and the potential strength 
more closely. At Vme = 0 , the peaks denoted by the blue symbols in Fig. 10.21 can only 
appear if they obey the selection rule An = 1 . After an external modulation potential is 
applied, peaks from other selection rules (An = 0 and An = 2 ) appear as shown by the 
black symbols. Under a weak modulated potential {Vme < 0.005 70 ), the peak frequency 
is hardly affected. If the electric potential continues to grow, some peaks start to split. 
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Figure 10.20. The optical absorption spectra for the composite field Bum = 5 T combined 
with Rmm = 395 and Vme = O.O 270 , where the contribntions from the transitions 

and kmid kmid are shown in (b) by the dotted and dashed curves respectively, (c) 
A comparison between the absorption spectra for Bum = 5 T combined with Rmm = 395 
and Vme = O.OO 570 and Bum = 5 T combined with Rmm = 1580 and Vme = O.O 270 . 
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e.g., Under a strong potential {Vme > 0.02 70), more peaks are 

created by other selection rules, such as An = 3 (square markers). Even for the transitions 
between band-edge states, peaks can be developed by the rule An 7^ 1 (solid triangles). 
Finally, the absorption frequencies related to the kmid states decrease rapidly with respect 
to the increment of Vme- This is due to the fact that the conduction and valence LSs at 
kmid move closer to the Fermi level. These theoretical predictions could be examined by 


the optical-absorption spectroscopy methods. 76,83-86 
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CO (To) 



Figure 10.21. The dependence of absorption frequencies, and on the 

modulated strength Vme- The absorption frequencies with the selection rules |An| = 
|n — n'l = 1 and |An| ^ 1 are indicated by the blue and black colors respectively. 
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8 Conclusion 


The results show that monolayer graphene exhibits the rich optical absorption spectra, an 
effect being controlled by the external ffelds. Such fields have a strong infiuence on the 
number, intensity, frequency and structure of absorption peaks. Moreover, there would exist 
the dissimilar selection rules for different external fields. In the presence of the uniform 
magnetic field, the magneto-optical excitations obey the specific selection rule An = 1, 
since the simple relationship exists between the two sublattices of a- and 6-atoms. As to a 
modulated magnetic field, an extra selection rule An = 0 is obtained, due to the complex 
overlapping behavior from two subenvelope functions in the wave function. However, the 
wave functions exhibit irregular behaviors under the modulated electric potential. As a 
result, it is difficult to single out a particular selection rule. In the composite ffelds, the 
symmetry of the LL wave functions is broken by the introduce of two kinds of modulation 
fields, the modulated magnetic and electric fields, a cause resulting in the altered selection 
rules in the absorption excitations. The one case is the uniform magnetic field combined 
with the modulated magnetic field. The extra selection rules, e.g., An = 2 and 3 , come 
to exist when Bmm is comparable to Bum- Another case is the uniform magnetic field 
combined with the modulated electric potential. The extra selection rules, e.g., An = 0, 2 
and 3 , would be generated in the increase of Vme- 

For the other graphene systems, the magneto-optical properties corresponding to a 
perpendicular uniform magnetic field deserve a closer investigation. For example, the AA- 


stacked bilayer graphene is predicted to exhibit two groups of absorption peaks; 72 however, 
the selection rule An = 1 is same as that of MG. As for the AB-stacked bilayer graphene, 
there exist four groups of absorption peaks and two extra selection rules (An = 0 and 2). 
The few-layer graphenes are expected to display more complex magneto-optical absorption 
spectra, mainly owing to the number of layers and the stacking configuration. 

In this chapter, the generalized tight-binding model is introduced to discuss mono- 
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layer graphene under five kinds of external fields. The Hamiltonian, which determines the 
magneto-electronic properties, is a giant Hermitian matrix for the experimental fields. It 
is transformed into a band-like matrix by rearranging the tight-binding functions; further¬ 
more, the characteristics of wave function distributions in the sublattices are used to reduce 
the numerical computation time. In the generalized tight-binding model, the 7r-electronic 
structure of MG is solved in the wide energy range of ±6 eV, a solution proving valid even 
if the magnetic, electric or composite field is applied. Moreover, the important interlayer 
atomic interactions, not just treated as the perturbations, could be simultaneous included 
in the calculations. The generalized model can also be extensible to other layer stacked 


systems, i.e., AA-, AB-, ABC-stacked FLGs 87-91 and bulk graphite. |57 ,92-941 
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Figure Captions 


Fig. 10.1 The primitive unit cell of monolayer graphene (a) in the absence and (b) in 
the presence of external fields. 

Fig. 10.2 (a) The Landau level spectrum for the uniform magnetic field Bum = 5 T. 
The Landau level wave functions corresponding to (b) the a- and (c) the 6-atoms. 

Fig. 10.3 The optical absorption spectra for (a) Bum = 5 T and (b) Bum = 10 T. (c) 
The dependence of the absorption frequency on the square root of field strength Bum- 

Fig. 10.4 The energy dispersions and the illustration of optical excitation channels for 
the modulated magnetic field along the armchair direction with Rmm = 500 and Bmm = 10 

T. 


Fig. 10.5 The wave functions of Quasi-Landau levels at (a)-(f) the original band-edge 
state ki with the quantum numbers rf = 1 and rf = 2, (g) and (h) the split point 

k 2 with rf = 1, and (i)-(l) two extra band-edge states k^ and k/^ with rf = 1. 

Fig. 10.6 The optical absorption spectra for (a) Bmm = 10 T at a fixed periodic length 
with modulation and polarization along the armchair and zigzag directions and (b) different 
modulation periods and field strengths with both the modulation and polarization along 
the armchair direction. 

Fig. 10.7 The dependence of the absorption frequency on (a) the period Rmm and (b) 
the square root of field strength Bmm- 

Fig. 10.8 The energy dispersions for the modulated electric potential along the armchair 
direction with Rmm = 500 and Vme = O.OSyo. 

Fig. 10.9 The optical absorption spectra corresponding to Fig. 10.8, which includes the 
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contributions from the ji and v states, respectively. 


Fig. 10.10 The optical absorption for (a) Vme = O.OSyo at a fixed periodic length with 
modulation and polarization along the armchair and zigzag directions and (b) different 
modulation periods and field strengths with both the modulation and polarization along 
the armchair direction. 

Fig. 10.11 The wave functions at different band-edge states, the and states for 

i = n. 

Fig. 10.12 The energy dispersions for (a) the uniform magnetic field Bum = 5 T by the 
red curves and the composite field Bum = 5 T combined with Rmm = 395 and Bmm = 1 
T by the black curves, (b) Bum = 5 T combined with Rmm = 395 and Bmm = 5 T, and 
(c) Bum = 5 T combined with Rmm = 395 and Bmm = 40 T. All modulated fields are 
applied along the armchair direction. 

Fig. 10.13 The wave functions with rf^'^ = 0 and rf = 1 oX for (a)-(d) the uniform 
magnetic field Bum = 5 T by the red curves and the composite field Bum = 5 T combined 
with Rmm — 395 and Bmm = 1 T by the black curves, (e)-(h) Bum = 5 T by the red 
curves and Bum = 5 T combined with Rmm = 395 and Bmm = 5 T by the black curves, 
and (i)-(l) Bum = 5 T combined with Rmm = 395 and Bmm = 40 T. 

Fig. 10.14 The optical absorption spectra corresponding to (a) the uniform magnetic 
field Bum = 5 T by the red curve and the composite field Bum = 5 T combined with 
Rmm = 395 and Bmm = 1 T by the black curve and (b) the composite field Bum = 5 T 
combined with Rmm = 395 and Bmm = 5 T. 

Fig. 10.15 The optical absorption spectra corresponding to the composite field Bum = 5 
T combined with Rmm = 395 and Bmm = 40 T by the red dotted curve and the pure 
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modulated magnetic field Rmm = 395 and Bmm = 40 T by the black curve. 

Fig. 10.16 The dependence of absorption frequencies, cu””' and cu””' with |An| = 
|n — n'l = 1, on the modulated strength Bmm- 

Fig. 10.17 The energy dispersions for (a) the uniform magnetic field Bjjm = 5 T 
by the red curves and the composite field Bum = 5 T combined with Rmm = 395 and 
Vme = O.OO 570 by the black curves, (b) Bum = 5 T combined with Rmm = 395 and 
Vme = 0 . 0270 . All modulated fields are applied along the armchair direction. 

Fig. 10.18 The wave functions for (a)-(d) = 0, = 1, rf = 5, and = 5 at 

kmid, (i)"(ii) = 0, = 1 and = 5 at and (l)-(m) = 5 at The results 

corresponding to the uniform magnetic field Bum = 5 T, the composite field Bum = 5 
T combined with Rmm = 395 and Vme = O.OO 570 , and Bum = 5 T combined with 
Rmm = 395 and Vme = O.O 270 are indicated by the red, black, and blue curves respectively. 

Fig. 10.19 The optical absorption spectra corresponding to (a) the uniform magnetic 
field Bum = 5 T by the red curve and the composite field Bum = 5 T combined with 
Rmm = 395 and Vme — O.OO 570 by the black curve, (b) the transitions of ^ k^^^ 
and kmid kmid for Bum = 5 T combined with Rmm = 395 and Vme — O.OO 570 by the 
dotted and dashed curves respectively, and (c) the higher frequency region of 0 . 1 ~ 0 . 27 o. 

Fig. 10.20 (a) The optical absorption spectra for the composite field Bum = 5 T 
combined with Rmm — 395 and Vme = O.O 270 , where the contributions from the transi¬ 
tions — 7 - and kmid kmid aro shown in (b) by the dotted and dashed curves 

respectively, (c) A comparison between the absorption spectra for Bum = 5 T combined 
with Rmm = 395 and Vme = O.OO 570 and Bum = 5 T combined with Rmm = 1580 and 
Vme = 0 . 0270 . 
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Fig. 10.21 The dependence of absorption frequencies, and on the 

modulated strength Vme- The absorption frequencies with the selection rules |An| = 
|n — n'l = 1 and |An| ^ 1 are indicated by the blue and black colors respectively. 
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